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ABSTRACT 



Radiances obtained from the 11 -y window channel by MRIR scanning 
radiometer of NIMBUS II were investigated in order to determine areas 
of time change of T 2 , which may be related to effective cloud top level 
of the five-day mean cloud state. A statistical model was developed 
that is useful in obtaining a valid specification for effective cloud 
top change, using a stepwise regression program. 

The parameters AT-j , AZ, and T^ appear to provide the necessary 
spatial and temporal stability to employ the specification scheme over 
a time span of ten days. These three parameters are all highly descrip- 
tive of AT 2 in terms of basic atmospheric dynamical relationships. It 
is recommended that a test be made of this hypothesis using a 24-hour 
specification procedure. 
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I. INTRODUCTION 



The diagnosis of weather features by sequential analysis of satellite 
maps of infrared radiances sensed during the lifetime of the NIMBUS II 
satellite, presents opportunities for understanding tropical and sub- 

I 

tropical atmospheres generally, and in particular the environment of the 
subtropical jet stream. These radiances were obtained by the MRIR 
scanning radiometer of NIMBUS II, which scanned in the 11-p window 
channel (channel 2) and gave an equivalent blackbody temperature (here- 
after called T^). This radiometer also scanned in the 6.7-y water vapor 
channel (channel 1), which corresponds to the effective water-vapor 
temperature (hereafter called T-j ) in the upper troposphere. For a more 
detailed discussion of MRIR data, the reader is referred to the NIMBUS 
il User’s Guide (Staff, 1966). 

Grid point maps of the two equivalent blackbody temperature fields 
T-| and T^ for the period 15 May through 28 July 1966 were provided in 
the form of five-day means through the kind auspices of the Laboratory 
for Atmospheric and Biological Sciences of NASA at Greenbelt, Maryland. 
The use of five-day mean maps presented consecutive fields with no data 
gaps on the quasi-globe defined between 70S-70N latitude (Fig. 1), 
although all of the transient daily features have undoubtedly been 
smoothed out. 

The purpose of this research was to specify statistically the time 
change of T 2 , denoted hereafter as AT^, especially in data-sparce areas 
over the oceans. The actual value of AT 2 may be related to the change 
in the effective or mean cloud-top level of the five-day mean cloud 
state at grid point (i,j). In turn the field of AT 2 (i,j) will be 
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statistically related to certain derived products which are descriptive 
of basic atmospheric processes. 



Areas of significant AT-] values are known to be related to areas of 
variation in the cirriform cloud-top level and/or variation in depth of 
moist layers associated with updrafts or downdrafts in the layer 6-10 km 
[Raschke and Bandeen (1967), Beran, et al . (1968)], During the course 
of this investigation it became evident that AT^ is a most descriptive 
parameter for specifying AT^, and hopefully is amenable to prognostica- 
tion for periods up to 12 hours in advance of the latest radiometric 
swath. If so, the method described here could be altered to provide a 
12-hour prediction for AT 2 - 

Although radiometric data covered the region 70S-70N in Mercator 
grid coordinates, the area of interest for this study was restricted to 
the limits 14. 3N to 56. 8N for which arse fivs-dav 5 BO mb contour dr.tr, 



were available. At the same time, these limits provided sufficient 
latitude to accomodate each of the subtropical jet streams considered. 
The longitude span of the subtropical jet area was restricted for the 
purposes of this study to the range 180W to 1201J. This area could 
normally be classified as data-sparce. During the data period of this 
study, the area displayed well defined five-day mean subtropical jet 
coverage. 
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Mercator grid point chart o/er which radiometric data were distributed. 



II. PROCESSING OF DATA 



The radiometric data were distributed over a Mercator grid point 
chart having i=l,...,73 grid points in the longitudinal direction, and 
j=l,...,41 rows of latitude covering the meridional span 70S-70N. The 
Mercator base chart had a superimposed square mesh grid, with a mesh 
interval ^equivalent to five degrees of longitude, true at the equator 
(Fig. 1). The radiometric data T(i,j) were subjected to an initial 
filtering which eliminated waves of two gridmesh lengths, but left 
unaffected all waves of length L j> 8cL 

As a possible dynamic predictor, the five-day mean 500 mb Northern 
Hemisphere non-overlapping contour fields for five-day map periods 
synoptic with the radiometric temperature fields were also obtained in 
puncned car a format covering the longitude dimension just described, 
but limited to 18 gridpoint rows poleward of latitude 14.3N (i.e., to 
56. 8N in Fig. 1). It should be noted that the contour map fields were 
sufficiently smooth that preliminary filtering was not considered neces- 
sary for these fields. 

Certain operations are defined as representing the time-mean and 
time-differencing operations on temperature data at each grid point of 
the successive maps. These time centered operations are defined below: 
Temperature-mean considered at time k: 

2T k (i,j) = T k + -|/2^”^ + T k-l/2^’^ 



Temperature-time difference at k.: 

AT k (i,j) = T k+ ]/2^»j) ‘ T k-l/2^’^ 
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where T^ = radiometric five-day mean temperature centered at time k+1/2 
i = x-direction grid indicator 
j = y-direction grid indicator 

k = 1,2,..., 14, the time-step indicator (half steps are 

centered at times of individual five-day means. Fig. 2). 



K*o x /z. 



k*l 



Fig. 2. Illustrating time-sequencing of five-day mean periods, and 
the time-index scheme used in centered time differencing. 



A similar procedure was followed to form AZ by subtraction of successive 
five-day mean 500-mb contour fields which are also centered at the half 
time-steps indicated in Fig. 2. 

A computer program was generated to carry out these time-centered 
compu Id Lions over the meridional range i4.3il-bo.8N as well as over the 
longitudinal extent of this investigation. This resulted in a set of 
fourteen consecutive maps, each of which consisted of five fields of 
time-means and time-differences of the following parameters: aT-j ^ , 

?! |^, AT^ ^ AZ^. These fields were programmed for production 
both in punched card and print-out forms. 

Fifteen consecutive subtropical jets in the eastern Pacific were 
manually determined using the fifteen original five-day mean 500-mb 
contour maps for the periods centered at times t=k+l/2, with k=0,...,14. 
The selection of individual jet-latitude coordinates was made by examining 
the test area 180W-1201J for the latitude of the tightest contour gradient 
that could be assigned nearest to the subtropical anticyclone on the 
five-day mean 500-mb chart. 
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These five-day mean subtropical jets lay well within the previously 
described latitude range of 14.3N-56.8N in which range most of the large 
contour-change AZ centers could be located. Furthermore the fifteen 
geostrophic mean jets varied only slightly in time and space, so that an 
average of six grid points lay both to the north and south sides of the 
jet. In order that the statistical test area be consistently more 
extensive than the jet area, the data set was expanded longitudinally to 
include two gridpoints west of 1801J and three to the east of 120W. Thus 
the sample size for each of the k=l,...,14 map-differenced times was 
constant at 216 in this investigation. 
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III. STATISTICAL PROCEDURES 



A. EXPERIMENTS FOR DETERMINING PREDICTION TIME-STEPS 

The Biomedical stepwise regression program, BMD02R [Dixon, 1966], 

was used in the statistical experiments described in the following pages. 

| 

A brief summary of symbols and equations follows. 

When the statistical specification is made for AT 2 , the predictors 
are understood to be X-j , X 2 , Xg, X^, defined at mean map time k, and each 
grid-point (i,j) as: 



X 1 ■ 


( T 1 ,k+l/2 " T 1 ,k-l/2^ 


= AT i.k 


X 2 - 


( z k+1/2 “ 1 k-1/2^ 


* AZ k 


X, = 


^2 k+1/2 + ^2 k-1/2^ 


* 2T ? 


X 4 

AT 2,k " 


^ T l, k+1/2 + T 1 ,k-l/2^ 
^ T 2 , k+1/2 " T 2 ,k-l/2^ 


' 2T l,k 



The variables defined in (1) appear in the sets of specification 
equations associated with AT 2 (k) in the form 

at 2 = a q + a 1 x 1 + a 2 x 2 + a 3 x 3 + a 4 x 4 (2) 

where A. are the regression coefficients, selected in a stepwise 
manner by the program BMD02R. 

In order to estimate the variable AT 2 centered at any one given 
map-time k from the concurrent predictor set in (1), it was necessary 
to determine how large a sample base of dependent data was needed. This 
question was tentatively answered by performing the following three 
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statistical experiments, each of which tested not only the significance 
of the specifications but also the stability of the regressions for at 
least one map-period. 

1. Each of thirteen time-centered periods (i,e,. at time k=l ,2, , ,13) 
was used to generate a specification equation (2) which could be used 
with the data of the next period (i.e., at times k=2,3, . . . ,14) leading 
to a set of thirteen verifications. This experiment gave the best indi- 
vidual period specifications and also single-period verifications on the 
test used in the experiments described below. The order of entry of the 
predictors varied less from case to case, implying greater stability in 
this system of deriving the specification equations than in tests 2 or 3 
described below. 

2. The seven odd-numbered map-time fields were pooled in order to gener- 

r) ^ipnlo rnmnncitp cnori f iVati on + ThG COCf f 1 Cl ^ c qt 

equation were used to operate on the variables X-| , X^, X^, and X^ of 
equation (2) at the even-numbered time periods, so that, in effect, a 
composite single-period forecast scheme was devised. In this case, the 
specification equation itself did not afford sufficient explained 
variance of the predictand aT^ for a valid comparison with the results 
of tests 1 and 3. The transgenerated equation explained even less 
significance when applied to the composited seven-day test data (at times 
k=2,4,. .. ,14) . This method indicates the lack of statistically homogene- 
ous samples when drawn from a meteorological period encompassing 13 
five-day mean maps, or 65 days. 

3. The final test undertaken was to use a pooled data base consisting of 
three consecutive periods (e.g., k=l,2,3) in order to generate a three-map 
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composite specification equation for use in predicting AT^ for the single 
following period (e.g., k=4). This set of apparent single-step forecasts 
was verified against the known data for time k=4. 

Similarly, the succeeding three-period pooled data base for times 
centered at k=2,3,4 was used for a composite statistical sample to pro- 
vide a verification at map-time k=5, and likewise on to the final period 
k=14. 



Test 3, which used a data base of pooled data from three consecutive 
maps in order to transgenerate an estimated AT^ for the fourth in 
succession actually involved on the average a full two-period verification 
test. This test was simply not as effective in describing explained 
variance on the independent sample as test number 1 described above, 
which involved single-map time-steps. Hence all results to be considered 
in Tables 2 and following a»'e ha«.<=>d upon the results cf test 1 only. 

As previously noted these comparative tests were primarily based upon 
the relative stability of the verification procedures carried out upon 
independent data. 

The definitions X-j , X^, X^ and X^ serve to simplify the description 
of the results in Table 2(a) below. This table refers to the order of 
entry of the four variables in the specification equation associated with 
AT£ of equation (2) . 

The determination of the F-statistic upon entry, computed at each 
step of the multiple regression program (BMD02R), forms the basis for the 
order of entry of the variables X^ in Table 2. 

The definition of the F-statistic at the Jth step of entry is defined 



as: 



Fj[l,n-J-1] = 



[cum. % exp. var.,step J] - [cum. % exp, var. , J-l] i 
[remaining % var. after step J] 
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Here is the number of grid points in the data sample. Since signifi- 
cance is sought at the 95% level of confidence for the Jth variable 
added, Fj is compared with a critical value Fj suggested by Miller (1962) 
as 

F J = F a/(P.-J+1 )^ ’ n " J+1 )’ a " ,05 

where is the number of predictors among employed in the 

r 

regression. The sequence of Fj values is easily deduced from tables 
of F-values at the indicated rejection limits for a/(P_-J+l). These 
critical Fj-values are listed below for the Jth selection step. 

TABLE 1. The critical Fj statistics at the 95% confidence limit 
for entry of the Jth predictor of a total predictor set P=4 
variables. The sample size is n=216 in each case. The last row 
shows the corresponding critical values when only three predictors 
are considered for entry by the stepwise reqression. 





J=1 


J=2 


J=3 


J=4 


P=4 


6.34 


5.83 


5.12 


3.89 


P=3 


5.83 


5.12 


3.89 





B. TEST FOR AN OPTIMAL NUMBER OF PREDICTORS 

The stepwise regression technique applied to AT 2 and its predictor-set 
X-| , . . . ,X^ on map-day 1 led to the results shown in Table 2(a). Comparison 
of the F , upon entry with the critical values in Table 1 at steps J=l,...,4 
indicates that the fourth predictor entering does so with a very low 
F-value and a resulting low probability of rejecting the null hypothesis 
that A^ equals 0 in equation (2). However the first three predictors all 
entered with confidence levels in Table 2(a) well in excess of the 95% 
level of belief. 
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Part (b) of Table 2 contains the specification results for AT^ at 
map-time 1 with the variable 2 T, excluded from the predictor set by 
means of an option of the BMD02R program. The predictor set in Table 2(b) 
is now reduced to AT-j , AZ and 2 T 2 . An objective rationale for the 
statistical preference of 21 ^ to 2T-| is given in more detail in IV(C), p 32. 

In Table 2(c), the stability of the specification equation at time k=l 
is tested by correlating the observed AT 2 values of sample map 2 against 
the values estimated from the transgenerated equation for AT 2 (k=2), made 
by using the day 1 coefficient matrix listed in Table 2(a). Symbolically 
this operation is performed using the estimator matrix equation: 



Here coefficients are determined form the program BMD02R and are 
listed in column 6 of Table 2(a). 

Some question exists regarding the relative stability of the four- 
predictor equation (4) compared to the three-predictor equation inherently 
given in Table 2(b). This three-predictor transgenerated estimator is 
denoted by the symbol AT 2 ( 2 ) for day 2 and has the form: 



Here the B^ are the coefficients listed alongside the predictors in 
Table 2(b). 




(4) 




(5) 
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TABLE 2. Specification of the predictors and their order of entry 
in equation (2) and for the map period corresponding to time k=l. 

Part (a) refers to the four predictor specification of AT 2 ; part (b), 
to the three predictor specification of AT^. Part (c) refers to the 
verification of A 12 ^ rorn equation (4) applied to test time k=2. 

Part (d) refers to the verification of aT^ from equation (5) applied 
to time k=2. 



1 Step 
J no. 


Variable 

entered 


Mult. Corr. 
Coeff. , 
step J 


Fj upon 
entry 


Coeff. of 
Jth var. 


% exp. 
var. 


J=1 


AT i 


0.7763 


324.625 


1 .20781 


60.27 


J=2 


AZ 


0.8149 


38.915 


0.02242 


66.41 


J=3 


2T i 


0.8322 


19.600 


-0.14271 


69.25 


J=4 


2T 2 


0.8328 


0.712 


-0.01269 


69.35 


o 

c 

V' 


ant Tprm - 






"»/i ' ’ # V LH 




J 








/ l * f / ^ v 






j=i 


at i 


0.7763 


324.625 


1 .22177 


60.27 


J=2 


AZ 


0.8149 


38.915 


0.02214 


66.41 


J=3 


2T 2 


0.8259 


12.028 


-0.04019 


68.21 


Constant Term - 






22.1258 




J=1 


A 

AT 2 


0.6543 


160.212 


0.93228 


42.81 


J=1 


AT 2 


0.6995 


205.054 


0.95789 


48.93 
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The verification of a prediction on an independent sample-set such 
as that of map-time 2, is judged to be significant when a large F-j- 
statistic results for either AT 2 or AT 2 compared to F-j , the critical 
value necessary at the 95% level (Table 1, with J=l). In other words, 
the question to be answered is whether the verifying F-statistic used 
upon independent test-data is greater than the critical F required for 
the 95% confidence level, i.e., is F(1 ,216) >_ F a ( 1 ,216) with probability 
1-a >^0.95? If the indicated F-statistic is sufficiently large, then 
there is a probability p >_ 95% of a linear relationship between A^ and 
the test variables (4) or (5). In parts (c) and (d) of Table 2, both 
F-statistics are sufficiently high to refute the probability of a chance 

✓N-/ 

relationship, although AT 2 (based upon equation (5)) is found to have a 
higher F-value than that of AT 2 of (4), 

C. REQUIREMENTS FOR A USEFUL SPECIFICATION FOR AT £ 

A useful specification will be defined as one satisfying the two 
criteria listed below. 

In the first place, the specification equation (4) or (5) at each 
step J must be considered significant at the 95% confidence level using 
the F-statistic upon specification: 

F(j ,n-J-l ) = % mean sq. var. exp! . by (4) or (5) 

% mean sq. var. unexpl . after use of (4) or (5) 

The F-statistic (J=1,...,P) of (6) is very nearly identical to the 
9 

computation of Rj, where Rj is the multiple correlation coefficient at 
the Jth step of entry in Table 2(a,b). The F-statistic of (6) is signifi- 
cant at a confidence level in excess of 95% whenever the F-statistic at 
J=P exceeds the fiducial limit F : 
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( 7 ) 



F(P,n-P-l) > F (P ,n-P-l ) , a =,05 
— a 

This requirement is conservatively satisfied at the final step J=P with 
n=216, whenever Rp 0.50. 

A AV 

The second criterion for the test variables aTq and/or AT 2 must 
also satisfy the inequality (7) with P=1 , for there now exists only one 
way of applying the specification equation to the verification test data 
of time k+1 . This latter criterion makes a requirement upon the time 
stability of the specification equation, a concept already used in III (A) 
in deciding upon the sampling procedures to be used. 

The second criterion may be most conveniently judged by computing 
the simple correlation coefficients R(AT 2> AT 2 )p +1 or R(AT 2 ,AT 2 ) k+1 , 
depending upon whether the four- or three-predictor regression (4) or (5) 
is being tested. The statistics are shown, for example in column 4 of 
Table 2 (c,a) as they pertain to map 2. 

D. ASSESSMENT OF THE USEFULNESS OF THREE-PREDICTOR RELATIVE 

TO FOUR-PREDICTOR CASES 

The test used in criterion 2 of II 1(B) was made at every verification 
time k=2,...,14 with the appropriate coefficient matrices of (4) and (5), 
respectively. It was found that the three-predictor equation led to the 
majority of stable verifications based upon consistently larger values 

A 

of the test-statistics R(aT 2 ,aT 2 ) compared to R(aT 2 ,aT 2 ). This is 
equivalent to saying that the three-predictor equation (5) gave a 
consistently higher F-level upon verification than did the four-predictor 
equation (4). For the reasons just cited and for economy of display, 
only the tabulations of specification and prediction of AT 2 (the three- 
predictor mode) will be given in IV, which follows. Furthermore, the 
useful rule for a satisfactory specification, R ^ 0.5 will be adhered to. 
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at least in 
significant 



the sense of providing a guide line to a statistically 
prediction. 
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IV. STATISTICAL RESULTS AND INTERPRETATIONS 



A. TABULAR LISTINGS OF THE INDIVIDUAL SPECIFICATIONS AND VERIFICATIONS 
The results of the specification of AT^ for map 1 and the application 
to the prediction of map 2 have already been listed in Table 2. Hence- 
forth the tabular listings v/ill be based upon the use of the three-predictor 
equations for both specification and verification. As an example, the 
results of the specification-prediction process for the case of map 2 
proceeding to map 3 are listed in Table 3. 



TABLE 3. (a) Specification of AT? by means of three predictors and 

their order of entry in equation (2) for the map corresponding to time 
k=2. (b) the predictand verification for aT ^ at time k=3, tested using 

fp3ncr!onov>afoH oninf i np ( vri+-h rhp rno-P-Pir lontc y'nc ij 1 +• -j nq -py«om TT|^ H 9 



Step 

no. 


Variable 

entered 


Mult. Corr. 
Coeff . , 
step 0 


Fj upon 
entry 


Coeff. of 
0th var. 
after 0=3 


% exp. 
var. 


0=1 


AT, 


0.6480 


154.920 


1 .25588 


41 .99 


0=2 


AZ 


0.7051 


32.748 


0.01502 


49.72 


0=3 


2T 2 


0.7140 


5.466 


-0.03471 


50.99 


Constant Term — 






18.67230 




0=1 




0.6049 


123.476 


1 .01990 


36.59 


Constant Term — 






0.54359 
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TABLE 4. (a) Specification of AT„ by means of three predictors and 

their order of entry in equation \2) for the, map corresponding to time 
k=3. (b) The predictand verification for AT 2 at time k=4, tested using the 

transgenerated equation (5) with the coefficients resulting from map 3. 



! Step 
no. 


Variable 

entered 


Mult. Corr. 
Coeff . , 
step J 


Fj upon 
entry 


Coeff. of 
Jth var. 
after J=3 


% exp. 
var. 


J=1 


AT i 


0.5846 


111.1034 


1 .30967 


34.17 


J=2 


AZ 


0.6514 


30.5582 


0.01701 


42.43 


J=3 


T 2 


0.6626 


5.5326 


0.03524 


43.90 


Constant Term -19.08235 




J=1 


A f 2 


0.6664 


170.9637 


1 .18358 


44.41 


Constant Term -0.11233 





TABLE 5. (a) Specification of AT 2 by means of three predictors and their 

order of entry in equation (2) for the, map corresponding to time k=4. 

(b) The predictand verification for Al '2 at time k=5, tested using the 
transgenerated equation (5) with the coefficients resulting from map 4. 



Step 
no . 


Variable 

entered 


Mult. Corr. 
Coeff. , 

Step J 


Fj upon 
entry 


Coeff. of 
Jth var. 
after J=3 


% exp. 
var. 


J=1 


AT i 


0.6423 


150.3109 


1 .63211 


41.26 


J=2 


AZ 


0.7037 


34.8372 


0.03241 


44.25 


J=3 


T 2 


0.7151 


7.0636 


-0.04216 


51.14 


Constant Term -- 






23.47249 




J=1 


iT 2 


0.5810 


109.0053 


0.83761 


33.76 


Constant Term -• 






-0.22636 
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TABLE 6. (a) Specification of AT 2 by means of three predictors and their 

order of entry in equation (2) for th§ map corresponding to time k=5. 

(b) The predictand verification for AT 2 at time k=6, tested using the 
transgenerated equation (5) with the coefficients resulting from map 5. 



Step 

no. 


Variable 

entered 


Mult. Corr. 
Coeff . , 
step J 


Fj upon 
entry 


Coeff. of 
Jth var. 
after J=3 


% exp. 
var. 


J=1 


fiT i 


0.4798 


64.0106 


1 .28656 


23.02 


J=2 


AZ 


0.6168 


51.6145 


0.04833 


38.04 


J=3 


T 2 


0.6175 


0.3010 


-0.00859 


38.13 


Constant Term -- 






4.17122 




J=1 


AT 2 


0.5753 


105.8482 


0.72407 


33.09 


Constant Term -- 






0.14901 





TABLE 7. (a) Specification of Al 2 by means or three predictors and their 

order of entry in equation (2) for the, map corresponding to time k=6. 

(b) The predictand verification for AT 2 at time k=7, tested using the 
transgenerated equation (5) with the coefficients resulting from map 6. 



Step 

no. 


Variable 

entered 


Mult. Corr. 
Coeff. , 

Step J 


Fj upon 
entry 


Coeff. of 
Jth Var. 
after J=3 


% exp. 
var. 


J=1 


4T i 


0.5615 


98.5142 


0.97251 


31.52 


J=2 


f 2 


0.6417 


34.9574 


-0.06165 


41.18 


J=3 


AZ 


0.6666 


12.4533 


0.01840 


44.44 


Constant Term -- 






33.53210 




J=1 


AT 2 


0.4619 


58.0505 


0.87025 


21.34 


Constant Term -• 






0.71712 
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TABLE 8. (a) Specification of AT 2 by means of three predictors and their 
order of entry in equation (2) for the, map corresponding to time k=7. 

(b) The predictand verification for AT 2 at time !<=8 > tested using the 
transgenerated equation (5) with the coefficients resulting from map 7. 



Step 

no. 


Variable 

entered 


Mult. Corr. 
Coeff . , 
step J 


Fj upon 
entry 


Coeff. of 
Jth var. 
after J=3 


% exp. 
var. 


J=1 


AT i 


0.4481 


53.7527 


0.95690 


20.08 


J=2 


AZ 


0.5919 


49.0300 


0.05084 


35.03 


J=3 


f 2 


0.6090 


6.9456 


0.04418 


37.09 


Constant Term 


— - 23.46455 




J=1 


A-/ 

iT 2 


0.6588 


164.0839 


1 .24020 


43.40 


Constant Term -- 






-2.38063 





pec rn cation or oy means or tnrce predictors 
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TABLt 9. (a) Sp 

order of entry in equation (2) for the map corresponding to time k=8. 
(b) The predictand verification for AT 2 at time k=9, tested using the 
transgenerated equation (5) with the coefficients resulting from map 8. 



Step 

no. 


Variable 

entered 


Mult. Corr. 
Coeff. 
step J 


Fj upon 
entry 


Coeff. of 
Jth var. 
after J=3 


% exp. 
var. 


J=1 


4T i 


0.5594 


97.4484 


1.26707 


31.29 


J=2 


AZ 


0.6643 


48.9321 


0.05761 


44.12 


J=3 


T 2 


0.6666 


1.2026 


0.02007 


44.44 


Constant Term -- 




-12.24481 




J=1 


AT 2 


0.5439 


89.9256 


0.43759 


29.59 


Constant Term -* 






0.75761 
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TABLE 10. (a) Specification of AT^ by means of three predictors and their 
order of entry in equation (2) for the map corresponding to time k=9. 

(b) The predictand verification for KT 2 at time k=10, tested using the 
transgenerated equation (5) with the coefficients resulting from map 9. 



Step 

no. 


Variable 

entered 


Mult. Corr. 
Coeff . , 
step J 


Fj upon 
entry 


Coeff. of 
Jth var. 
after J=3 


7o exp. 
var. 


J=1 


AZ 


0.3859 


37.4493 


0.02426 


14.89 


J=2 


AT i 


0.5368 


41.6514 


0.69409 


28.81 


J=3 


T 2 


0.5785 


14.8342 


0.06375 


33.47 


Constant Term -- 


34.89899 




J=1 


2f 2 


0.1449 


4.5912 


-0.17341 


2.10 


Constant Term -- 






1 .31401 





11 ( Z* \ i r A-P AT ^ K w m a in r q T -f- Ia v->/A /A r\ v*/a rl a 'f' a *a rJ f Ko 1 ’ ^ 

order of entry in equation (2) for the map corresponding to time k=10. 

(b) The predicted verification for AT 2 at time k=ll, tested using the 
transgenerated equation (5) vrith the coefficients resulting from map 10. 



Step 
no . 


Variable 

entered 


Mult. Corr. 
Coeff. , 
step J 


Fj upon 
entry 


Coeff. of 
Jth var. 
after J=3 


% exp. 
var. 


J=1 


T 2 


0.3431 


28.5494 


-0.09331 


11.77 


J=2 


AT 1 


0.3480 


0.8302 


0.06180 

. j 


12.11 

L . 


J=3 


AZ 


0.3506 


0.4407 


0.00469 


12.30 


Constant Term -■ 






52.30643 




J=1 


AT^ 


0.1015 


2.2278 


-0.10086 


1.03 


Constant Term - 






1 .17698 
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TABLE 12. (a) Specification of AT 2 by means of three predictors and their 
order of entry in equation (2) for the map corresponding to time k=ll. 

(b) The predictand verification for AT 2 at time k=12, tested using the 
transgenerated equation (5) with the coefficients resulting from map 11. 



Step 

no. 


Variable 

entered 


Mult. Corr. 
Coeff . , 
step J 


Fj upon 
entry 


Coeff. of 
Jth var. 
after 0=3 


% exp. 
var. 


J=1 


AT, 


0.4562 


56.2325 


0.96817 


20.81 


J=2 


T 2 


0.5057 


13.6216 


-0.08226 


25.57 


J=3 


AZ 


0.5298 


7.3782 


0.02377 


28.07 


Constant Term -- 




43.16838 




J=1 


AT 2 


0.5698 


102.8863 


1.14851 


32.47 


Constant Term 4.59990 j 



TABLE 13. (a) Specification of AT 2 by means of three predictors and their 
order of entry in equation (2) tor the map corresponding to Lime k~i2. 

(b) The predictand verification for ST 2 at time k=13, tested using the 
transgenerated equation (5) with the coefficients resulting from map 12. 



Step 

no. 


Vari able 
entered 


Mult. Corr. 
Coeff. , 
step J 


Fj upon 
entry 


Coeff. of 
Jth var. 
after J=3 


% exp. 
var. 


J=1 


AT, 


0.6457 


153.0436 


1.57667 


41 .70 


J=2 


AZ 


0.7031 


32.6091 


0.02136 


49.44 


J=3 


T 2 


0.7114 


5.0147 


0.04039 


50.61 


Constant Term -20.36539 




J=1 


AT ^ 


0.6688 


173.2133 


0.76579 


44.73 


Constant Term -1.61974 
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TABLE 14. (a) Specification of AT 2 by means of three predictors and their 
order of entry in equation (2) for the_ map corresponding to time k=13. 

(b) The predictand verification for AT? at time k=14, tested using the 
transgenerated equation (5) with the coefficients resulting from map 13. 



Step 

no. 


Variable 

entered 


Mult. Corr. 
Coeff . , 
step J 


Fj upon 
entry 


Coeff. of 
Jth var. 
after J=3 


% exp. 
var. 


J=1 


AT i 


0.5876 


112.8828 


1.00762 


34.53 


J=2 


AZ 


0.7372 


92.4117 


0.04857 


54.34 


J=3 


T 2 


0.7399 


1.9182 


0.02137 


54.75 


Constant Term -- 






-11 .92155 




J=1 


at 2 


0.6030 


122.2797 


1 .01096 


36.36 


Constant Term -- 






1.40262 





TABLE 15. (a) Specification of AT? by means of three predictors and their 
order cf entry in equation (2) for the, map corresponding Lo time k“14. 

(b) The predictand verification for 2T? a t time k=l , tested using the 
transgenerated equation (5) with the coefficients resulting from map 14. 



Step 

no. 


Variabl e 
entered 


Mult. Corr. 
Coeff. , 
step J 


Fj upon 
entry 


Coeff. of 
Jth var. 
after J=3 


% exp. 
var. 


J=1 


AT i 


0.5100 


75.2309 


1.12036 


26.01 


J=2 


AZ 


0.5982 


32.4248 


0.03983 


35.79 


J=3 


T 2 


0.6485 


22.9468 


0.09734 


42.06 


Constant Term -- 






-52.64392 




J=1 


ar 2 


0.7011 


206.8751 


0.72137 


49.15 


Constai 


it Term -- 






-0.57721 
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In addition it is desirable to list all corresponding sets of aT£ 
specification and verification for each of the twelve subsequent pairs 
of days. For example, Table 4 is analogous to Table 3, except that the 
initial specification is for time k=3 and the verification using the 
transgenerated equation is then applied to k=4. The remaining Tables 
5-15 are listed in sequence in view of a similarity in the pattern of 
the statistics presented in the listings, when useful specification- 
predictions occur. In contrast, a distinct change in this pattern occurs 
when a specification-prediction fails, an example of which occurs at 
map-time k=10. 

In order to exploit the data for maximum usefulness, the specifica- 
tion from k=14 was used as a predictor for k=l , the results of which are 
shown in Table 15. 

B. A bUMMARY OF THE MEAN SPECIFICA1 ION-PREDICTION kESULTS 

A statistical parameter that clearly summarizes the specification- 
verification results for the two-map statistical sequences is the final 
percentage explained variance. Hence it is advantageous to list this 
statistic for each specification map together with its corresponding 

verification. This is done in Table 16. 

~2 

The parameter R is identified with the mean overall explained 
variances in the set of twelve useful specifications and corresponding 
verifications, which appear at the end of columns 2 and 7. Then R, the 
mean coefficient of determination, has the value R" = 0.6911 and .6082 at 
the indicated stages, respectively. 

Table 16 also lists in column 3, for each specification map k, the 
percentage unexplained variance, which is given by (1-R ). 
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TABLE 16, Summary of percentage explained variances of AT ? 
at the specification step k, and at step (k+1 ) in the 
verification process. 



Map- 

time 

k 


% Exp. 
Var. AT 2 
on map k 


% Unexp. 
Var. for 
map k 


Standard 

Deviation 


Standard 

Error 


Map- 

time 

k+1 


% Exp. Var.' 
AT 2 on 1 
map k+1 


1 


68.21 


31 ,79 


5.30951 


2.9736 


2 


48.93 


2 


50.59 


49.41 


5.20436 


3.6400 


3 


36.59 


3 


43.90 


56.10 


5.43297 


3.9325 


4 


44.41 


4 


51 .14 


48.86 


6.40811 


4.4229 


5 


33.76 


5 


38.13 


61 .87 


5.82407 


4.5978 


6 


33.09 


6 


44.44 


55.56 


4.83820 


3.5816 


7 


21.34 


7 


37,09 


62,91 


5.RQ^nQ 


a. 6610 


8 

1 


43.40 

1 ' » 


8 


| 44.44 


55.56 


6.62934 


4.3862 


9 


29.59 




9 


33.47 


66,53 


5.75243 


4.1150 


10 


2.10 


10 


12.30 


87.70 


5.14620 


4.8410 


11 


10.92 


11 


28.07 


71.93 


5.15380 


4.2230 


12 


i 

32.47 


12 


50.16 


49.84 


5.65208 


4.0115 


13 


44.73 


13 


54.75 


45.25 


5.12623 


3.4371 


14 


36.36 


14 


42.06 


57.94 


6.32972 


4.8106 


1 


49.15 


Mean : 


45.11 

i 


54.89 


5.6710 


4.2171 




36.99 
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In addition, the table lists the standard deviation a of AT 2 
and its standard error (a^) f° r eac ^ time k=l,...,14, The star 
errors in each case may be computed from equation (8): 




standard error 

standard deviation 



= 1-R 



(8) 



It should be noted that the two "unreliable" prediction cases 
(10 and 11), which occur because of a single anomalous specification 
at k=10, have not been included for averaging purposes in Table 16. 



C. NORMAL PATTERNS OF PREDICTOR-ENTRY FOR USEFUL SPECIFICATIONS OF AT 2 

The statistical feature most apparent among the useful (R ^.0.5) 
specification cases listed in Tables 2 through 15 was the order of entry 
of the predictors. When the multiple correlation coefficient had values 
in excess of entry of O.b, the order ot entry in these tables was 
consistently AT^ , aZ, and T^. As mentioned in 1 1 1 A , it is the value 
of the F-statistic upon entry which determines the order of entry, and 
in ten of the twelve cases considered to be useful, the predictors 
entered in this particular order. The other two cases of useful speci- 
fication corresponded to the data samples for maps 6 and 11, in both of 
which the order of entry was AT^ , T 2 , and AZ. For both of these cases 
the explained variance for specification was acceptable due to the large 
F upon entry at the first step, i.e., that which introduced the variable 
AT 1 . 

An "unreliable" assessment was given to the specification equation 
for map 9 (Table 10) for which the irregular order of entry (AZ, aT-j , T 2 ) 
and associated relatively small F upon entry values led to coefficients 
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which could not be meshed with those entering the specification equation 
for map 10, Thus map 10 gave poor verification of map 9, and also poor 
specification when introduced for the verification of map 11. The 
instability of these two map-sets may be seen at once by reviewing 
Tables 10 and 11, with particular attention given to the coefficient 
matrix of the variables AT-j , AZ, T^. Note particularly the systematic 
change in the matrices shown relative to those entered immediately before 
and after maps 9 and 10. 

Using percentage explained variance as a guide to the verification 
of maps 10 and 11, it is apparent from Table 16 that the map-sequence 
10-11 should be deleted as useful specification-prediction combinations. 
When this has been done, the mean fractional explained variance upon 
specification and verification are obtained: 0.4511 and 0.3699, 
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al explained variance at both specification and verification, namely 0.25. 

At this point a brief discussion of the option of selecting rather 
than T"i and is applicable. The primary reason for this decision was 
the high value of the linear correlation coefficient, r^.T-j), in all 
useful specifications. As an example, the correlation matrix based upon 
map time 2 is shown in Table 17 below. 



TABLE 17. The predictand-predictor 
linear correlation matrix for map 2. 





AT 2 


ATt 


AZ 


l2_ 


Tl 


AT 2 


1 .000 


.648 


.124 


.141 


.041 


AT 1 


.648 


1.000 


-.227 


.363 


.230 


AZ 


.124 


-.227 


1.000 


1 

O 

00 


-.266 


T 2 


.141 


.363 


-.048 


1 .000 


.768 


A 


.041 


.230 


-.286 


.768 


1.000 
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Note in particular that F(T| ) = 0.768 so that there exists 
considerable redundancy of similar information when both and are 
retained in a four-variable specification. The decision to eliminate 
T-| rather than by the control-delete option of BMD02R is then based 
upon the small standard deviation of the observed T , , relative to that 
of T 2 (approximately in the ratio of 1 to 3). Furthermore, using the 
CONDEL option in testing the three predictor equation, which included 
T| rather than the deleted T 2 , resulted in a slight shrinkage of the 
multiple correlation coefficients upon verification. 
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v. EXTENSION OF THE RESULTS TO 24-HOUR SPECIFICATIONS 



A. RE-EXAMINATION OF THE FIVE-DAY SPECIFICATION 
PROCEDURE FOR POSSIBLE EXTENSION 

In the preceding discussion, AT^ values centered at the five-day 
mean map time k were found to be linearly related to the independent 
variables AT-j , AZ and T 2 by means of a statistical equation of the 
following form: 

AT 2 = B q + B,AT 1 + B 2 AZ + B 3 T 2 (9) 

where each of the variables of the right side is likewise time-centered 
at map time k. In addition, the order of importance of the variables on 
the right side of (9), in thirteen of the fourteen cases, was significant 
using the step-wise regression of the BMD02R program. The ord^ was 
consistently AT-j , AZ, 

In order to further establish the stability of (9), a verification 
test was conducted, based upon the results at time k, to the five-day 
mean data at time k+1 . This resulted in only a slight shrinkage in 
percentage explained variance using the same coefficient matrix 
(Bq» > ^ 2 , B^) on the independent data of map k+1. 

This stability test was desirable to assess the usefulness of the 
diagnostic equation (9) in deriving a valid specification. In the 
process of testing (9) for statistical stability, unreliable specifica- 
tion in one case, k=10, led to two cases of inadequate verification, at 
times k— 1 0 and 11. Statistically this has been explained in IV as due 
to lack of homogeneity in the five-day mean coefficients B Q , B-j , B 2 , Bg, 
which occurred at map-time k=10, relative to the longer-term sample 
means of these coefficients over the full test-period of 14 maps. 
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A pertinent conclusion to be drawn, then, is that the statistical 
model evolved in IV is useful in obtaining a valid specification for 
the effective cloud-top level change for a five-day mean cloud state at 
grid points (i,j), using the Standard Atmosphere lapse rate of 6.5K(km)~^ 
Thus a negative change in ATp of 6.5K is equivalent to an effective cloud 
top ascent of 1 kilometer, if fully overcast conditions exist. 

On the other hand, if the mean cloud cover at an oceanic grid point 
remains constant at (T, where C < 1 . 0 , and an observed decrease in T 2 
occurs, then the ascent in the effective cloud tops (which covers only 
the fractional area C) is deductible from 

AT. 

AZj = 3 , y = 6.5 K/km 

K Y-C 



This formula follows since the window channel radiometer senses the 
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the effective cloud top level over the fractional area C of cloud-cover. 

Finally, the parameters AT^ , AZ and T 2 appear to provide the necessary 
spatial and temporal stability to employ the specification scheme over 
time spans extending as long as ten days in most of the cases listed 
(Tables 2,. ..,15). 



B. 24-HOUR STATISTICAL SPECIFICATION OF AT. AND THE PHYSICAL 
INTERPRETATION OF THE PREDICTOR VARIABLES 

The problem to be considered now is that of the analogous use of 
strictly 24 hourly data. It would be possible to use the BMD02R program 
to obtain a diagnostic equation of the form: 

AT 2 (24) = C Q + C^ + C 2 AZ + C 3 T 2 00) 

based entirely upon independent variables observed at 24 hourly intervals, 
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and time-centered twelve hours before the latest data report has been 
received. The assumption is that good daily mapped versions of the 
and T.j fields have been received in suitably calibrated digitized Tg B 
form from the designated processor in Washington (NESS), and communicated 
instantly to the satellite data-user in the field. 

The tests applied in this paper do not include 24-hour best-fit 
specification because AT^, AT^ , AZ and T^ were not available in synop- 
tically time-centered form at appropriate times. However, if useful 
specifications were derived on the basis of five-day averaged data of 
the same form, it is expected that the specification equations of at 
least the same or better efficiency should result for 24-hour averaged 
data. This concept will be discussed from the point of view that AT^ 
depends upon basic atmospheric physical parameters , which are reflected 
thvTMjnh AT a 7 pprl T d c (jc* H r ' ci 1 1 cT r • Moreover The Tirnr— 



sequencing of the physical relationships, and therefore of the specifi- 
cations, should be more precise for shorter periods (24 hours) than for 
longer period specifications, such as five-day periods. During this 
longer period, atmospheric diffusion tends to smear the effects of the 
physical relationships which bear on ATg- 

Equation (10) is therefore proposed as a significant regression to 
be applied to 24-hour mean data fields much as such regressions v/ere 
found to be descriptive for five-day means. Further supporting argu- 
ments below are for a 24-hour specification equation of form (10) based 
upon the physical usefulness of the variables on the right side of (10). 
1. The 24-Hour Time Chanae, AT, 

* — i 1 

Raschke and Bandeen (1967) have shown that properly calibrated 
values of T-j (when used in conjunction with T^) yield the mean relative 
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humidity in the tropospheric layer extending between six to ten kilo- 
meters. AT-j may therefore be interpreted as either a change in the mean 
relative humidity of this column or as a vertical velocity producing the 
same humidity effect. 

Beran, et al . (1968) showed that 24-hour positive changes in AT-| 
are associated with areas of convergence near the 200 mb jet level, with 
resultant downward forced vertical motion giving reduced mean humidities 
in the upper convergent tropospheric layer. The use of five-day means 
in the present investigation, as contrasted with that of Beran, resulted 
in the loss of some discrimination regarding AT-j , since the five-day mean 
fields of AT-j , near the mean polar front jet had much smaller gradients 
than were found by Beran et al . on a single-day change basis. 

2 • The 24-Hour Hejoht Change at 500 mb, AZ 

24-^O'jr h 0 *j n n f rhanno f '! o 1 He a + mK pmM*lr\\/or| l-iov'p cuiro 

they have been found to be useful in analysis of cloud-development seen 
in satellite pictures. In that context, Z is usually presented in its 
scale-pattern decomposition: 

7 = 7 + z 

L l sr l sd 

( 11 ) 

AZ =AZ sr+ AZ sd 

Here, since the long-wave field Z<^ is temporally conservative, it 
follows that AZ = AZ<~p , and therefore, 

AZ = Z $D (t=0) - Z $D (t=-24) . 

Distributions of AZ<.q are known to manifest the short wave 500-mb 
vorticity patterns which are superimposed upon the large-scale steering 
patterns, the latter being the Z<^ fields at 500 mb. In turn, the AZ^ 
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patterns are related to photographic features of cloud analysis, e.g., 
the comma-shaped cloud with maxima of positive vorticity advection. 

Thus AZ patterns occurring over the preceding 24 hours represent changes 
in the relative vorticity patterns over the area and are instrumental 
not only in effecting changes in the photographic depiction of cloud 
masses or bands but also changes in their vertical extent, reflected by 
parameter T 2 . In summary, the parameter AZ, as defined for the past 24 
hours, should be expressive of a cloud-structure change in AT^ resulting 
from the vorticity-advection change field within the grid area. 

3 . The Large-Scale Steering Parameter, Tp 

Finally, in scale decomposition applications, the 500-mb SR 
field is known to represent the steering component of the 500-mb wind 
field insofar as vorticity changes are concerned. Space-average contour 

hoinji'f' *f’i p"l rlc oni 11 VP 1 op ^ O "f’hp 7 . lwpv^p no'r a %/a *i 1 an *J p 

for any of the tests conducted in this study. Nevertheless, a time- 
smoothed 500-mb map is very nearly equivalent in its long-wave steering 
characteristic to the space-smoothed Z S p-maps which may be derived from 
an individual 500-mb map, at the time-centered date. 

The Northern Hemisphere five-day mean contour heights have been 
added to the data file for use in a BMD02R program for the same set of 
15 data periods as were originally investigated. This was done in a 
manner exactly analogous to that used in forming fields of T^, which is 
a predictand entering the right side of (9). The five-day Northern 
Hemisphere fields of T 2 , for the 1966 period of NIMBUS II, have been 
intensively studied by Professor F. L. Martin, who has found that they 
afford a striking resemblance to the corresponding five-day mean 500-mb 
contour fields, period-for-period north of 15° latitude. 
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This fact suggested the supplementary test of correlating the 
five-day means of T 2 and T over the test area of this investigation 
(216 grid points), based upon a pooled series of five consecutive maps 
of five-day means of both T 2 and 1. The ensuing regression analysis gave 
the result: 

Z(gpm) = 1919.0293 + 6.9477T 2 

with la linear correlation: R(Z, T 2 ) = .8506. 

This result gives a rationale for employing T 2 in the correlation 
originally tested in the specification equation (9) where, except for AZ, 
only radiometric variables were entered into the regression equation. 

It thus becomes clear that the five-day mean T 2 field has nearly the same 
smoothed properties as the five-day mean Z-field. 

It is proposed that, for the 24-hour specification, the use of 
the time-centered field of T 2 , when space-smoothed to grid points on an 
NMC synoptic-scale chart, reflects primarily the SR-advecting component 
of the 500 mb field. Thus in (10), the space-smoothed T 2 field should 
by analogy perform the role of the 500-mb SR field in steering the 
individual disturbances. For this reason the space-smoothed T 2 (from 
the 12-hour earlier swath) has been included as predictor in the 24-hour 
specification equation (10) for AT 2> 

Finally, it has been concluded, using both the references and 
tests cited here, including related knowledge from satellite-picture 
interpretations, that the three parameters AT^ , AZ, and T 2 which appear 
in (10), are likely to be descriptive of AT 2 - It is recommended that a 
test be made of this hypothesis using the 24-hour specification procedure. 

From the preceding discussion, each of the parameters just cited 
is representative of a large-scale dynamical factor affecting cloud-structure 
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in the atmosphere. While there is undoubtedly non-linear interaction 
between the physical parameters proposed, the use of a linear regression 
model tends to separate nonlinear effects into linearized ones. Hence 
the following summary of the statistical parameters included in (10) is 
presented here with the suggested physical interpretation which they 
represent: 

(1) aT-j is indicative of vertical motion change over the past 24 hours. 

(2) aZ is indicative of the relative vorticity propagation by steering 
and of development during the steering. 

(3) The space-smoothed field is indicative of the space-smoothed 
500-mb field. It is also indicative of a 24-hour mean radiometric 
estimate of the cloud-top level over oceanic grid points. In this 
later connection, T 2 ( i , j ) should also be related to average surplus 

r\ -f pf i \m n ^ n \ / i n moan ro] rimn n 4* ( n ^ pnn fho v>o _ 
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fore to the eddy-generation effect on an existing or incipient 
disturbance (Katayama, 1967). 
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VI. PROPOSALS FOR FUTURE STUDY ; A SUMMARY 



In order for an equation of form (10) to have predictive value for 
a period of 12 hours over a given test area, it will be necessary to have 
a polar-orbiting satellite such as NIMBUS IV, reporting T-j and T£ data 
in equivalent blackbody form at 12 hourly intervals. In this way, the 
24-hour interval centered 12-hours before receipt of the latest report 
would provide the resulting specification equation, and could then be 
applied to a time period centered now at time t=0, concurrent with 
receipt of the latest radiometric maps provided aT^ can be computed. 

This scheme would require the t=0 field of T£, smoothed to synoptic 
scale, as well as prognosticated fields of aT-j and aZ for the succeeding 
twelve hours. A 24-hour AZ prognosis is presently available from the 
Primitive liquation model. Prognostic fields of Ai-j are not presently 
available; however it is not unreasonable to suggest that with further 
tests and evaluation using satellite-provided water '-vapor channel fields 
such a prognostic scheme could be developed. 

NIMBUS IV included a temperature-humidity infrared radiometer (THIR) 
experiment, the purpose of which was to measure by channels 1 and 2 
daytime and nighttime surface and cloud-top temperatures and cloud 
mapping; by channel 1, atmospheric water vapor mapping (NIMBUS IV User's 
Guide, 1970). THIR data are archived in separate daytime and nighttime 
swaths. The problem is then one of obtaining prognosticated 12-hour 
fields of T-j , so that T^ (12) - T^ (-12) = AT is known for use in (10). 

If the temperature distribution is initially known, the vertical 
distribution as well as total amount of water vapor in a vertical column 
above 500 mb may be inferred (Raschke and Bandeen, 1967). The inference 
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of water-vapor absorber mass below 500 mb has not been very successful 
for radiances in the 6.7y channel because the effective-intensity emission 
in channel 1 arises from a level in the upper troposphere (Shenk and 
Salomonson, 1970), Furthermore, characteristic T-j configurations provide 
information about the upper-tropospheric wind field. Martin and 
Salomonson (1970) have used radiometric information from NIMBUS II MRIR 
to locate the core of the subtropical jet and deduce maximum wind speed. 

Thus since is indicative of in-cloud temperature and is subject to 
change with changing cloud-structure much as a satellite photograph 
indicates, the most favorable radiometric field to consider for prognosti- 
cation is that of T-j . A successful prediction in the T-j field is 
reasonable to hypothesize since T-j is descriptive of the mean water- 
vapor mass above 500 mb; that is, at a level generally above the top of 



+ tr»*i r] r! 1 O — Y' 1 O t ! H C 1»M ^ T-h ■? o r> v*r\ /-* ^ /■> 4- -i 4- r\ A T i H 



-P A r> 1 A r* +' A ! 



would become available for use with the known prognostic field of AZ, and 
the current field. Thus aT^, the most variable radiometric parameter 
of equation (10), would be amenable to computation by the regression 
method proposed here. 

This avenue of approach seems to offer a fruitful use of MRIR radiances 
as predictors of effective cloud-top change. Furthermore, the result 
would afford a three-dimensional aspect to cloud-distribution forecasting. 
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